Field-doping of Ceo crystals: Polarization and Stark splitting 
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We investigate the possibility of doping Ceo crystals by applying a strong electric field. For an 
accurate description of a Ceo field-effect device we introduce a multipole expansion of the field, the 
response of the Cgo molecules, and the Stark splitting of the molecular levels. The relevant response 
coefficients and splittings are calculated ab initio for several high symmetry orientations. Using a 
group theoretic analysis we extend these results to arbitrary orientations of the Ceo molecules with 
respect to the external field. We find that, surprisingly, for the highest occupied (HOMO) and the 
lowest unoccupied molecular orbital (LUMO), respectively, the two leading multipole components 
lift the degeneracy of the molecular level in the same way. Moreover the relative signs of the 
splittings turn out to be such that the splittings add up when the external field induces charge 
into the respective level. That means that when charge carriers are put into a level, its electronic 
structure is strongly modified. Therefore, in general, in Ceo field-effect devices charge is not simply 
put into otherwise unchanged bands, so already because of this their physics should be quite different 
from that of the alkali-doped fullerenes. 

PACS numbers: 73.61.Wp,73.90.+f,74.70.Wz 



I. MOTIVATION 

The proposal of doping Cgo crystals in a field-effect de- 
vice (FET) and the possibility of metallic conduction and 
even superconductivity in such devices had raised wide- 
spread interest. While the revelation of dishonest data 
handling in some casesi led to a severe damping of the ini- 
tial enthusiasm, fundamental aspects of field effect dop- 
ing remains a timely and interesting problem. For field- 
effect transistors made with self-assembled monolayers 
this question was addressed in Ref. and for the reported 
enhancement of the superconducting transition tempera- 
ture in Cpo crystals intercalated with haloform molecules 
in Refs.^J Attempts to observe the field-effect in graphite 
were reported in Ref.0. Here we address, from a theoret- 
ical point of view, the question how strongly Cgo can be 
doped in an electric field before its electronic structure is 
substantially changed and how this structure changes in 
even stronger fields. This is relevant for understanding 
the fundamental features of field-effect devices based on 
Cgo and involves a number of interesting physical prob- 
lems. 

It appears that doping Ceo crystals in a field-effect de- 
vice (FET) is very hard to achieve in practice, one of 
the reasons being the exceptionally strong fields that are 
required. Very strong fields, however, not only induce 
charge carriers, but also polarize the molecules, and, due 
to the Stark effect, in general, lift degeneracies. These 
effects are of particular importance, as Ceo is quite polar- 
izable, and as its molecular levels are highly degenerate. 
The term "field-doping" naively implies that these effects 
are small, such that the main effect of the external field 
is inducing charge carriers into electronic levels which are 
essentially unaffected by the field. It is clear that if the 
external field is strong enough, the electronic structure 
of the crystal will be strongly modified by the field, and 



one thus can no longer speak of doping. A fundamental 
question about field-effect devices is therefore connected 
with the doping levels achievable, before the electronic 
structure of the active material is substantially altered. 
Moreover we discuss in detail how the electronic structure 
is changed, when the external field is beyond this "dop- 
ing limit" . In this regime the commonly used analogy of 
the field-doped Cgo with the alkali-doped fullerenes no 
longer holds. Nevertheless a field-effect device with high 
carrier concentration would be an interesting device in 
its own right. One might, e.g., speculate that for Cgo in 
a strong electric field the electron-phonon coupling is en- 
hanced compared to unperturbed Cgo, as for molecular 
orbitals of lower symmetry less couplings are forbidden 
by symmetry 

For a first estimate of the fields involved in field- 
doping, consider a simple capacitor. Given a charge 
per area a on the plates, the electric field between the 
plates is E = Ana. Assuming that in field-doped Cgo 
the induced charge resides in the top-most layer^ in- 
ducing n elementary charges per molecules requires an 
external field (originating from the gate electrode) of 
E cxt = 27rn/A mo i, where A mo \ is the area per molecule 
in the top-most layer of the crystal. For a Cgo crys- 
tal with lattice constant a ~ 14 A typical areas per 
molecule are ^4(m) = V3a 2 /A for the (Iff)-plane and 
^(ooi) = a 2 /2 for the (OOl)-plane. Even though these ar- 
eas are quite large, the external fields necessary for field 
doping are substantial, being of the order of I V/A per 
induced elementary charge per molecule. This is, how- 
ever, not the field experienced by a molecule. As Cgo 
is highly polarizable (a ~ 83 A 3 ), in the solid, the ex- 
ternal field at the site of a molecule is screened by the 
polarization of the neighboring molecules: A monolayer 
of dipoles p centered on the lattice sites Ri, generates 
a field Edip = — pY^i^ol-R-i ~ i?o| — 3 at Ro, where the 
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sum is over all sites in the monolayer, except Rq. For 
the (111) layer, this sum is about 31.2/a 3 , for the less 
dense (001) layer it is about 25.6/a 3 . As the dipole mo- 
ments p of the molecules are induced by the screened field 
E SCI = E cxt + Edi p at the site of the molecule (p — aE SCT ), 
we find, by self-consistent solution, that the external field 
is reduced by about a factor of two. 

Inducing charge and polarizing the molecules is, how- 
ever, not the only effect of the external field. In addition 
it also leads to a splitting of the molecular levels — the 
Stark effect. As the molecular orbitals of Ceo have a def- 
inite parity, a homogeneous field splits the levels only in 
second order. Thus for low fields the splitting is quite 
small, but increases quickly for larger fields. We can ex- 
pect that the splitting of the molecular levels disrupts the 
electronic structure of the crystal only when it becomes 
of the order of the band width, which is about 1/2 eV 
in Ceo- Calculations indicate that the splittings of the 
T\ u - and iJ«-levels in a homogeneous field are surpris- 
ingly small, being less than 1/2 eV up to E SCT w 2 V/A 
(cf. figure EJ). Given the crudeness of the argument, it 
thus seems entirely feasible that doping of a few elemen- 
tary charges per molecule could be achieved before the 
electronic structure of the Ceo is substantially changed. 

In the argument above we have described the C6o 
molecules as polarizable points. It is then natural to 
refine the model by considering also higher multipoles. 
Such a multipole expansion is particularly suitable for 
Cgo as the molecules are nearly spherical. The approach 
is then the following: First we calculate the response of 
a Cgo molecule to external multipole fields. Then we 
self-consistently solve the electrostatic problem for a lat- 
tice of molecules in a homogeneous external field. This 
provides us with the multipole expansion of the screened 
field acting on each molecule in the solid. Given that 
screened field, we can then determine the splitting of the 
molecular levels as a function of the induced charge. 

The organization of the paper reflects this approach: 
In section [H] we describe the density functional calcu- 
lations for determining the multipole response and the 
Stark splitting for a Ceo molecule in a multipole field for 
several symmetrical configuration. Using group theory, 
in section IIIII the irreducible parameters for the multi- 
pole response are determined. This allows the calculation 
of the polarization for arbitrary configurations. In sec- 
tion IIVI we give an analogous treatment for the Stark 
splittings and explicitly show how the splitting changes 
as the molecule is rotated relative to the external field. 
In section [V] we use these ingredients to self-consistently 
solve for the screened field seen by a molecule in a charged 
monolayer. The splitting of the molecular levels in this 
self-consistent multipole field and the effect of this split- 
ting on the density of states is presented in section IVll 
Our conclusions are presented in section IVlTI The meth- 
ods for calculating the response and splitting for an ar- 
bitrarily oriented external field from the results of the 
density functional calculations that were performed only 
for special orientations are described in the appendices. 




FIG. 1: Orientations of the Ceo molecule. The coordi- 
nate system is chosen such that the x-axis is pointing to 
the right, the j/-axis upward, and the z-axis towards the 
reader. Thus, from left to right, the molecules are oriented 
with their 2-, 3-, and 5-fold axis along z. The different ori- 
entations are obtained by rotating the molecule about the 
y-axis: by arctan(2 — r) to go from the first to the second, 
and by arctan(r) to go from the first to the third orientation. 
t — (V5+ l)/2 is the golden ratio. 

Appendix ^ gives an example of how to calculate the 
response of a molecule using the irreducible parameters 
derived in section ITTT1 In appendix 151 we derive the cou- 
pling matrices needed for the calculating the level split- 
ting when the molecule is rotated in the external field. 
Finally, appendix |0 gives the derivation of the matrix 
describing the field generated by a lattice of identical 
multipoles at the origin, which is needed for finding the 
self-consistent electrostatic field. 



II. RESPONSE OF A C 60 MOLECULE 

To determine the response of a Ceo molecule to exter- 
nal multipole fields we have performed all-electron den- 
sity functional calculations using Gaussian-orbitals^The 
basis set comprises 5s4p3d for carboni and we use the 
Perdew-Burkc-Ernzcrhof functional. 8 

In our calculations we apply an external multipole field 
and study the change in the multipole moments of the 
charge density and the splitting of the molecular levels 
as a function of the strength of the external field. To take 
advantage of the molecule's symmetry we consider mul- 
tipole fields with the z-axis along the 2-, 3-, and 5-fold 
axis of the molecule (cf. figure 0. As these axes are each 
contained in a mirror plane, which we chose to be the x-z- 
plane, we can treat the thus oriented molecule as having 
symmetry group C^v, C^ v , and G§ v respectively. Apply- 
ing external multipole fields with I > 0, this symmetry is 
maintained if the fields are proportional to the real part 
of the spherical harmonic Y\ m , where m is an integer mul- 
tiple of the order [n =2, 3, or 5) of the symmetry axis. 
Likewise, the response of the charge density will only have 
multipole components proportional to Re(Y; m ), with m 
an integer multiple of n. Calculations were done for such 
symmetry conserving multipole fields up to I = 6. For 
the 3-fold axis oriented along z, we have, in addition, 
calculated the response to external fields proportional to 
Im(y"5 3), Im(Y6,3) 5 and Im(l6 i.e., with a symmetry 
lowered to C3 . As we are interested in the linear response, 
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FIG. 2: Splitting of the H u and T\ u levels of a C60-molecule 
in an external multipole field V(r) = Vi m 7?* m (r) with z along 
a 2-fold axis. Due to parity, for small external fields, the split- 
ting for (Zm) = (10) is quadratic in the external field, while 
for (7m) = (20) it is linear. The V; m are given in atomic units 
and R g — 7 bohr. (Im) = (10) corresponds to a homoge- 
neous field with E z = — Vio, thus Vio Rg = 1 in atomic units 
corresponds to a homogeneous field of about 1 V/A. 



geometries can be extended by group theory to arbitrary 
orientations of the molecule relative to the external mul- 
tipole field. 

All calculations have been performed for the equilib- 
rium geometry of the unperturbed, neutral Ceo molecule. 
To estimate the effect of an external field on the shape of 
the molecule, we have relaxed the structure in the pres- 
ence of homogeneous external fields of up to 1 V/A. We 
find only small changes (up to about 0.005 A) in the 
lengths of the bonds (1.40 A for the short and 1.45 A for 
the long bonds). Likewise, the polarizability changes by 
less than 1.5%. 

Finally, we have calculated the total energy of the iso- 
lated Ceo ion as a function of its charge q (spin unpolar- 
ized calculation with relaxed geometries) and extracted 
the second order term 1/2 U$q 2 . To compare with pre- 
vious calculations, 9 we find for the polarizability (mul- 
tipole field with I = 1) of the neutral molecule a =9.3 
bohr 3 /atom, and Uq =3.2 eV. 



III. IRREDUCIBLE RESPONSE PARAMETERS 

The polarizability a of a molecule describes the linear 
dependence of the induced dipole moment p = a E on 
the applied electric field E. For a multipole expansion, a 
becomes a matrix ai imi i 2 „ l2 describing the response to 
all multipole fields^ To fix the notation (which follows 
Ref. ITo|) we briefly review the definition of the multipole 
response matrix. 

The solutions of the Laplace equation V 2 V(r) — are 
given by 



we have considered small multipole fields and made sure 
that the calculated response of the charge density is in- 
deed proportional to the strength of the external field. 
We find that the linear response of a Cgo molecule is 
very similar to that of a metallic sphere of radius 4.4 A. 
This effective sphere radius shows a slight increase with I. 
In addition there are weak off-diagonal terms. To judge 
the accuracy of our calculation, we have checked how well 
the selection rules, that are not already imposed by the 
C nv symmetry, are fulfilled for these off-diagonal terms. 
From the calculated response, we have determined the ir- 
reducible linear response coefficients, which will be given 
in the next section. 

In addition to polarizing the molecule, the external 
field also splits the degenerate molecular levels of the Ceo 
molecule. As the unperturbed molecular orbitals have a 
definite parity, for a multipole field with odd I there will 
be no first-order splitting — the quadratic Stark effect. 
On the other hand, a multipole field with even I can cou- 
ple states of like parity, so in that case the splitting is 
linear. This is shown in figure |2 In section IIVI it will 
be demonstrated how the splitting of the HOMO- and 
LUMO-orbitals that were calculated for high-symmetry 



V (r) = V c + V 1 = ]T Vi m Rt m (r) + Q lm V lm (r) , (1) 

1 tn 

where the two terms denote the external potential (V°), 
and the induced potential (V 1 ) due to a charge distribu- 
tion p located around r = 0. Note that both, the Laplace 
equation as well as the expansion of V^r) into multipoles 
only holds for r which lie outside the charge distribution. 
We have introduced the regular and irregular spherical 
harmonics 



Special cases for the regular spherical harmonics are 
i?oo = 1 an d Rio = z, hence, the external field, Vbo = Vq 
corresponds to a constant shift, and Vio = — E z is the 
z-component of the electric field. For the irregular spher- 
ical harmonics we have loo = 1/r and Iio = z/r 3 , thus 
for the induced potential, Qoo = 9 gives the monopole 
charge while Qio = p z is the dipole moment. Generally, 
the coefficients Qi m are the multipole moments of the 
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charge distribution p 

Qim = J d 3 rp(r)R lm (r). (4) 

Decomposing the charge distribution p = p + Ap into 
the unperturbed charge density and the change in the 
charge density due to the external potential, we obtain a 
decomposition of the multipole moments Qi m = Qi m + 
AQi m . Within linear response, the coefficients AQi m of 
the induced multipole moments depend linearly on the 
coefficients Vi m of the external potential, which defines 
the linear-response matrix ai imi i 2 m 2 - 

(5) 

i2"l2 

where the sign takes into account that the induced fields 
oppose the external fields. Then a\ mi im 2 gives the dipo- 
lar response tensor, while aoooo is the self-capacitance 
U . 

The interaction energy of the molecule with the exter- 
nal potential is E cxt = J d 3 r p(r) V c (r), which, using the 
previous definitions, reduces to 

E ext = J2v lm Qt n . (6) 

im 

Therefore, Vi m and Qi m are pairs of conjugate variables 
and the total energy of the molecule as a function of the 
external field is given by 

Etot = —~ m imi l 2 m 2 Vl imi VJ* m2 + ^2 Vim (Qlm)*- 

l 2 m 2 

(7) 

Since this is a quadratic form, we see that the matrix 
a is hermitian. We can make a real and symmetric by 
unitarily transforming to a real basis (using ^/2R,e(Yi m ) 
and ^/2lm(Yi m ) instead of Y/ m and Yi.-m for m ^ 0). 

The structure of the response matrix a depends on 
the symmetry of the molecule. For a metallic sphere 
of radius R the response is isotropic, i.e., a is diago- 
nal in the basis of the spherical harmonics: oti imi i 2 m 2 = 
Si ± i 2 5 mim2 R 2ll+1 . Lowering the symmetry to icosahe- 
dral, J/j, introduces some anisotropy. To understand the 
response matrix for Ceo we have to consider how the irre- 
ducible representations (IR) of the rotation group 50(3) 
split into IRs of the Ih (see table QJ. An external mul- 
tipole field of angular momentum I can only give rise 
to a response with angular momentum I', if both IRs 
of the 50(3) contain a common IR of the 1^. In par- 
ticular, because of parity, fields with even (odd) I can 
only give rise to responses with even (odd) V . Further- 
more, as the irreducible representations with I < 2 are 
also irreducible with respect to the for / < 2 we have 
a imim' — a iS mm '- Thus, restricting the multipole ex- 
pansion to / < 2, the response of Cqo is isotropic, with 
a Q fn 8.1bohr, an w 556bohr 3 , and a 2 « 44100 bohr 5 . 
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TABLE I: Decomposition of the irreducible representations 
(IR) of the rotation group 50(3) into the IR of the icosahedral 
group I h C 50(3). 



For I > 2 the space spanned by the spherical harmonics 
Yi m is no longer irreducible with respect to the Ih- Thus 
we need to find linear combinations of the spherical har- 
monics that span the irreducible representations of the 
icosahedral group. We call them Yi x k where I and x de- 
note the IR of the 50(3) and Ih, respectively, while the 
index k labels the functions within an irreducible rep- 
resentation of the Ih- If in the decomposition an IR x 
should occur several times, we would have to introduce 
an additional multiplicity label. However, as can be seen 
from table up to I = 6 each IR appears at most once. 
We therefore suppress the multiplicity label here. Ex- 
plicit expressions for the basis functions Yi x k, can, e.g., 
be found in Ref. IT2I chapter 16, or Ref. Il3l table 4.2. 

In the new basis, the matrix a is built of blocks of 
diagonal matrices 

ajixikj l 2 x 2 k 2 — aijzixi) S XlX2 Sk ± k 2 , (8) 

where ai 1 i 2 (xi) constitute the minimal set of parameters, 
and, a being real symmetric, a; 1 ; 2 (a;i) = a^i^Xi). The 
matrix elements of a were calculated up to I = 6 using the 
results of the density functional calculations described in 
section ITfl The a^^ixi) are listed in table |nj From 
this minimal set of independent parameters we can de- 
termine the response for arbitrary orientations of the Ceo 
molecule relative to the external multipole field. An ex- 
ample of how to do this is given in appendix^ The prac- 
tical advantage of this procedure is clear: We only need 
to perform density functional calculations for a number of 
highly symmetric configurations, for which the numerical 
effort is much reduced. Using group theory the response 
for arbitrary configurations can then be determined from 
these special cases. 

Note that the group theoretical approach presented in 
this section is particularly elegant in the case of a neutral 
molecule which has icosahedral symmetry. Upon charg- 
ing, orbitals become partially filled and the symmetry is 
reduced which leads to a higher number of irreducible 
response coefficients. Furthermore, the symmetry of a 
charged molecule depends on how the additional charge 
arranges in the degenerate orbitals. This is a subtle ques- 
tion in the case of an isolated Cgo molecule and involves 
Jahn- Teller effects and Coulomb interaction in competi- 
tionii. In the present work we restrict the analysis to 
neutral molecules. 
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aoo(A s ) 


1.019 
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aii(Ti u ) 
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ai5(Tiu) 
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022(^9) 
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a 6fi (H g ) 
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Q44(-n 9 J 
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„ (ft \ 
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0.074(6) 


a 2 e(H g ) 


-0.023(13) 






Ot46(H s ) 
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a 4 6(G 9 ) 


-0.337(13) 



TABLE II: Linear response coefficients of a neutral Ceo 
molecule derived from the results of our density functional 
calculations. The matrix elements ai 1 i 2 (x) are normalized 
with Ro — 8.25 bohr. By comparing matrix elements deter- 
mined from the response to potentials applied along the 2-, 
3-, and 5-fold axis, we have determined, wherever possible, 
the uncertainties in the values of the matrix elements. The 
value of aoo(^4 9 ) is given by the quadratic term of the change 
of the ground state energy upon charging of the molecule. 
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FIG. 3: Energy levels of the C60 molecule as calculated by 
density functional theory. The Fermi energy for the undoped 
molecule is indicated. 



IV. LEVEL SPLITTING 

In this section, a minimal set of parameters describing 
the Stark effect in a neutral molecule is deduced using the 
Wigner-Eckart theorem for the icosahedral symmetry^. 
This is achieved in the framework of density-functional 
perturbation theorjii& where an external perturbation V 
induces a change in the effective potential (self-consistent 
field) V which couples the degenerate energy levels. In 
linear order the change is given by 



AU ett (r) = V(r)+AK'(r) 



dn 



n—n(r) 



An(r), (9) 



where AV(r) = e 2 J dr'An(r')|r — r'| _1 is the change in 
the induced potential due to the linear change An(r) in 
the electron distribution. The last term is the exchange- 
correlation potential. Within linear response, the change 
in the charge distribution An as well as the induced po- 
tential AU 1 have the same symmetry (with respect to 
Ih) as the external perturbation V c . Furthermore, if 



n(r) is the unperturbed charge distribution, then the fac- 
tor 22s£| n _ , r \ in the last term of © has A g symmetry 
and does not change the symmetry of An. Consequently, 
AV has the same symmetry as V c and, using the no- 
tation of the previous section, can be written as 



AV 



clIV 



£ 

Ix k 



V lxk f l xk (r) + 0[V? x 



leh 



(10) 



where f l xk {v) are partner functions of the IR x of Ih- Note 
that f k contains spherical harmonics Yi> xk with all V al- 
lowed by table [I] The coupling of the degenerate levels 
is given by the matrix elements of AV with respect to 
the eigenstates of the unperturbed molecule which are 
denoted by \nxk), where n is the quantum number dif- 
ferentiating between orbitals with the same IR x. In this 
context, the functions f xk play the role of tensor opera- 
tors of Ih and the Wigner-Eckart theorem can be used to 
write the matrix elements as 



(n 2 x 2 k 2 \AV eS \n 1 x 1 k 1 ) 



(11) 



Ixk 



Vi x k ^2 *A(ni»i n 2 x 2 ;lx) Cl 2kl (A; x 2 x\; x). 

A 



The coefficients C k2ki (A; x 2 xi;x) denote the 3jm symbols 
(or Clebsch-Gordan coefficients) of Ih and are entirely 
determined by the icosahedral symmetry. In the present 
work they were taken from Ref. Il2| and are discussed 
in detail in appendix 151 In order for C k2ki (X;x 2 xi;x) 
to be non-zero, the IR x needs to be present in the de- 
composition of x\ ® x 2 . If a; occurs more than once in 
x\ <8) x 2 then the multiplicity index A is required. From 
this selection rule follows again that even potentials cou- 
ple linearly whereas odd potentials couple only in sec- 
ond order (quadratic Stark effect). Finally, the factors 
t\(n\X\n 2 x 2 \lx) are the coupling constants which con- 
stitute the minimal set of parameters describing the level 
splitting. 

We analyze in more detail the level splitting of the 
LUMO (x = T lu ) and HOMO (x = H u ) for I = 1 and 
I = 2 external potentials which correspond to an electric 
field and a quadrupole potential. We will see in Sect ion Ivl 
that these two multipoles are dominant in the case of a 
charged Ceo layer exposed to an electric field. Within the 
icosahedral symmetry Ih, the two potentials form part- 
ners for the IR T ltl and H g respectively (see tabic . 
Details of the calculation are given in appendix The 
coupling constant derived from the level splitting calcu- 
lated by DFT are given in table lllll In the case of the odd 
I = 1 potential, only second-order coupling to closest-by 
orbitals was considered (T l9 for the LUMO, H g and G g 
for the HOMO, see figure |3J, which, however, as shown 
below, gives very satisfactory results. For the splitting 
of the HOMO under the I — 2 (x — H g ) potential, there 
are two coupling constants because H g occurs twice in 
the product H u <g) H u = A g 8 T lg T 2g ® 2G g 2H g . 
In appendix IbI the coupling matrix H' T {H' H ) describ- 
ing the level splitting of the LUMO (HOMO), "due to an 
applied (Im) = (10) and (Im) = (20) potential along an 
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X2 
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t\(xi X2\ x) 




Tx 9 


T\ u 




-0.663 




Ti u 


Hg 




-0.225 


H u 


H g 


Tiu 




-0.730 




G g 


Tiu 




-0.730 




H u 


H a 


1 


-0.520 




H*i 


H a 


2 


-0.018 



TABLE III: Coupling constants t\(xi X2',l) as denned by 
equation 111! (indices m, n 2 , I are dropped) among the en- 
ergy levels near the Fermi energy. The coupling constants are 
given in units of eV (7bohr) !+1 e _1 where I = 1 for x = Ti u 
and I — 2 for x — H g . The negative sign of the coefficients 
is due to the negative charge of the electrons. Note that the 
coupling constants t(H u H g ;Ti u ) and t(H u G g ;Ti u ) are equal 
up to the third digit which is due to fact that the H g and G g 
levels have almost the same radial dependence because of the 
their closeness in energy (see Fig|3j. 
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arbitrary direction of the molecule, is calculated using 
perturbation theory. To a very good approximation, the 
result can be cast into the form 

H kA 9 ) = (KcL + Vxc^Cn^O), (12) 
H' Hu (0) = (v* di + V 20 da) C Hn (0) , (13) 

where Ci, C2, d\ and di are constants depending on the 
coupling constants of table II I II and the energies given 
in figure |21 Ct 1u (0) and Ch u {&) are matrices within 
the LUMO and HOMO subspace respectively and de- 
pend on the angle 6 between the z-direction of the 
(Im) = (10), (20) potentials and the 5-fold axis of the 
molecule, (for more details see appendix [Bj . In figure 0] 
and figure the splitting of the LUMO and HOMO is 
shown using the previous relations along with the points 
calculated by DFT. The group-theoretical fit is very sat- 
isfactory. Note that the splitting of the LUMO is in- 
dependent of the orientation of the molecule when only 
coupling among the LUMO or to the T\ g is considered 
(see appendix lB|l . Furthermore, relations (|12fl and 
are quite remarkable, as they imply that the contribu- 
tions of the I = 1 and I = 2 potential lift the degeneracy 
of the molecular levels in the same way, and thus the to- 
tal splitting is given by the sum of the I = 1 and I = 2 
splittings. 

V. MOLECULES IN A LAYER 

In this section a simple model of a Ceo FET is con- 
sidered. As mentioned above a FET can be understood 
as capacitance where one plate is the gate and the other 
plate is the material (here Ceo) which is investigated. An 
analysis of this device, in particular the calculation of the 
charge distribution, was done in a previous work^. It was 
found, that the charge concentrates on the first layer in 
the high doping regime. In what follows, we take this as 
a motivation to consider a single layer of Cgo molecules 



FIG. 4: Splitting of the LUMO for an (Zm) = (10) and 
(Zm) = (20) potential as the molecule is rotated by n/2 
about the y-axis. The orientations shown in figure are in- 
dicated. The splitting calculated with DFT at these high 
symmetry orientations are indicated by the filled symbols. 
The lines give the result of perturbation theory, fitted to the 
calculation for the 3-fold axis. Upper panel: (Zm) = (10) 
potential with Vio = 0.143 e (7bohr)~ 2 . The lower level is 
twofold degenerate. Lower panel: (Zm) = (20) potential with 
V 20 = 0.180 e(7bohr)" 3 . The upper level is twofold degener- 
ate. 



which acts as a plate of a capacitance. We use the re- 
sponse of a neutral molecule (table |nj to describe the 
electrostatic behavior of the molecules in the layer. The 
doped charge is taken care of by adding a monopole on 
every site. The molecules are then exposed to the electric 
field arising from the gate as well as the monopole fields 
of the neighboring molecules. In order to simplify the cal- 
culation, we assume a perfect lattice, either square, for 
the (001), or triangular, for the (111) plane of the fee lat- 
tice formed by the Ceo molecules in the bulk. In order for 
the sites to be equivalent we neglect the non-spherically 
symmetric part of the response of the Cgo-molecule and 
use the averaged response given by 

spher.-sym. r c \ * a hxk hxk (TA\ 

a l 1 m 1 l 2 m 2 — °hh °m 1 m 2 / , of | 1 ' 

x k 

where the sum is taken over all components of the l\ 
subspace. Since the response of Cgo for multipoles I < 2 
is isotropic (cf . section llll[l , this averaging is exact for 
I = 1,2. 

Because of translational invariance, the total potential 
is given by 

Utot(r) = U e (r)+^U i (r-Ri) ) (15) 

R.; 

where the sum is taken over all lattice sites R^. At a 
given site, say = 0, the total potential can also be de- 
composed as Vtot(r) = U scr (r) + U 1 (r) where the screened 



7 





-5.83 


> 


-5.832 






>> 


-5.834 


i— 

CD 








LU 


-5.836 




-5.838 




-5.82 


> 




SO, 


-5.84 


>, 




O) 




1— 
CD 


-5.86 


C 




LU 


-5.88 




-5.9 







~- — < 

1 i-l — i 1 1 1— 


» 


• 


— i 1 1 


C2v C3v , 


C5^^^ 
* C2v 



FIG. 5: Splitting of the HOMO for an (lm) = (10) and 
(lm)=(20) potential. See figure 0] for more details. In this 
plot the similarity of the splitting resulting from the different 
potentials (cf. eqn. is particularly striking. 



potential V scr (r) = V c (r) + £r^o V(r - R*) contains 
the external potential V c (r) as well as the sum of all 
induced potentials of the other molecules. This sum de- 
pends linearly on the induced potential V 1 and, hence, 
the coefficients of V SCI are given by 

+ E & 

imi l%mi Ql 2 m 2 - (16) 

l 2 m 2 

The matrix (3 is entirely given by geometry and discussed 
in appendixim On the other hand, the screened potential 
ysct j nc | uces a potential AV 1 as given in equation 
Therefore the total induced potential V 1 — V 1 ' + AV 1 is 





square 


triangular 


Im 






V lm" 




i n 
1U 


0.862 


0.499 


0.996 


0.530 


20 


-0.230 


-0.189 


-0.280 


-0.219 


30 





0.039 





0.053 


40 


0.0133 


0.0030 


0.0177 


0.0015 


4c4 


0.0174 


0.0097 







TABLE IV: Components of the bare and screened potential 
for the square and triangular lattice with a/y/2 — 10 A. The 
coefficients are in units of q/(7bohr) !+1 , where q is the charge 
per Ceo molecule. Note that 7bohr ~ 3.7 A is about the 
radius of the Ceo molecule. 




^ Square Lattice: 

V 6 "(0,0,z) 

V 6,, ([1-z 2 ] ,,2 ,0,z) 

Triangular Lattice: 

V e,, (0,0,z) 

V 6,, ([1-z 2 ] 1/2 ,0,z) 

0.5 

z/7a n 



FIG. 6: Screened potential V SCI across the layer with the 
parameters from Tab If VI for the square and triangular lattice. 
The gate is assumed to be on the left. The lower lines cor- 
respond to a cut through the center of the molecule whereas 
the upper lines are along half-circles with the radius of the 
molecule (7bohr « 3.7 A). 



Qhmi — Qlxmi 7 , C-limx ! 2 m 2 ^; 2 m2- (1^) 

Equations l|16|) and l|17(l can be combined by eliminating 
the coefficients Qi m which yields 

^ ^ E^+M^fam^X (18) 

i2"l2 

^irai = "T" y ' fthmi l 2 m 2 Ql 2 m 2 ' (1^) 

l 2 m 2 

where the VJ^"° describes the bare potential arising from 
the external potential (the electric field of the gate) 
and intrinsic moments of the molecules (the induced 
charge, i.e. monopoles). The square (triangular) lat- 
tice in the presence of the electric field has the rota- 
tional symmetry C± v (Cq v ). As a consequence, only com- 
ponents with these symmetries are non-zero and there- 
fore they are given by Ke(Yi rn ) with m a multiple of 4 



(6). Using relation i|18|) and i|19|) the screened poten- 
tial can be calculated. The non-zero components en- 
tering (fTT?f> are the monopole charge Qqu = 9 an d e l ec ~ 
trie field Vio = — Eoatc- As the FET is overall neutral 
^Gatc = — 27rq/A mo i with A(noi) = a 2 /2 for the square 
lattice and ^4(m) = ^/ia 2 /A for the triangular lattice. 
The results are given in tabic IfVI and graphically de- 
picted in figure The components (Im) — (10), (20) 
are most dominant and higher ones are at least one or- 
der of magnitude smaller. This justifies a posteriori the 
assumption of spherical symmetry because l|14|l is exact 
for (Im) = (10), (20). From figure[S]it can be seen that 
the electric field is efficiently screened within the layer. 
Note that decrease of the electric field yields negative 
sign of i>2o in table IIVI We also have checked the influ- 
ence of adjacent layers of Ceo and of a close-by dielectric 
(with a dielectric constant e — 10). The effects on the 
parameters in table ITVl were less than 2 %. The reason 
is that the field inhomogeneities induced by a 2D lattice 
of multipoles decay exponentially outside the lattice. 
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VI. SPLITTING IN SELF-CONSISTENT 
MULTIPOLE FIELD 

We are now in the position to estimate the effect of 
the external field on the electronic structure of the Ceo 
molecules in the monolayer that carries charge. To do 
so we have performed density functional calculations for 
a molecule in the self-consistent multipole fields as de- 
termined in the previous section (cf. table llVfl . Figure 
shows the splitting of the molecular levels in the self- 
consistent field for a (001) monolayer, where the molecule 
is oriented with one of its two-fold axes pointing in the 
direction of the external field. The maximum energy dif- 
ference between split states is given in table [V] and com- 
pared to the result from the perturbative formula l|12(l 
and i|13|) . which is in good agreement for \q\ < 2. As 
figure demonstrates, similar results are obtained for 
other geometries. Qualitatively, the results are also in 
agreement with an approximative tight-binding calcula- 
tion published earlier in Ref. ^3 As expected we find 
that the stronger the external field, i.e., the larger the 
induced charge, the stronger the splitting. We notice 
also a pronounced asymmetry in the splitting: when the 
monolayer is charged with electrons the splitting is dif- 
ferent from when it is charged with holes. Again, the 
reason is parity: Because of parity an external homoge- 
neous field, and more generally any multipole potential 
with I odd, leads to a quadratic Stark effect. Hence, 
for odd / the splitting is independent of the sign of the 
field. For I even, however, the levels split already in first 
order (linear Stark effect), so the splitting changes sign 
with the external field. Thus the asymmetry originates 
from the multipoles with even I. Moreover, because of 
the first-order versus second-order effect, even though the 
largest even multipole (20) is significantly smaller than 
the largest odd multipole (the screened external field) , it 
contributes considerably to the splitting. 

We have seen in section IIVI that for the HOMO and 
LUMO of C 60 the splittings caused by (10) and (20) 
multipole potentials are essentially additive (cf. equa- 
tion dj). Thus we expect them to add up for one sign 
of the external field, while they should compensate for 
the opposite sign. In fact, for electron doping the split- 
tings of the HOMO seem to almost perfectly cancel, while 
upon hole doping the splitting of the HOMO is essentially 
doubled compared to the splitting caused by the screened 
homogeneous field alone. Thus for the HOMO, the split- 
tings happen to add up for an external field that induces 
charge carriers into that orbital — a situation that is 
particularly unfavorable for hole doping. For the LUMO 
the situation is similar: For an external field that induces 
electrons in the LUMO the splitting is enhanced. So it 
turns out that the contributions of the higher multipoles 
conspire to enhance the splitting of the orbital that car- 
ries the induced charge. For both, HOMO and LUMO, 
the splitting becomes comparable to the respective band 
width for a field that corresponds to about two charge 
carriers per molecule. 





LUMO 


HOMO 


q 


DFT 


Pert. 


DFT 


Pert. 


-2e 


0.312 


0.305 


0.049 


0.001 


-le 


0.105 


0.102 


0.062 


0.060 

















le 


0.006 


0.002 


0.178 


0.179 


2e 


0.057 


0.097 


0.455 


0.476 


3e 


0.151 


0.297 


0.779 


0.892 



TABLE V: Maximum energy difference (in eV) between split 
HOMO and LUMO states, respectively, as a function of dop- 
ing q of the square lattice. 2nd and 4th column are the DFT 
results from figure 3rd and 5th column are calculated by 
perturbation theory as described in section IIVI 



The calculations reported above have been performed 
for an uncharged molecule. Considering instead the 
splitting for a molecule that carries the proper induced 
charge, the splittings are substantially reduced. This is 
easily understood: Due to the splitting, the electrons will 
fill only the lowest of the split levels. But since the in- 
teraction between electrons in the same orbital is larger 
than the interaction of electrons in different orbitals, the 
occupied levels will be shifted upwards, compared to the 
levels that were left empty - thereby reducing the split- 
ting. We are, however, not interested in the splitting per 
se, but in the effect of the splitting on the band struc- 
ture of a monolayer. Thus allowing electrons only in the 
energetically lowest levels would mean that in the lat- 
tice, the electrons are not allowed to hop to energetically 
higher ones of the split levels. This implies that the orig- 
inal band structure would already be separated into a 
set of bands originating from occupied, and another set 
of bands originating from the empty orbitals. To elimi- 
nate the undesired differences in the interaction between 
electrons in the split orbitals, we therefore work with the 
splittings obtained for a neutral molecule. 

To estimate the effect of the Stark splitting on the den- 
sity of states (DOS), we have performed tight-binding 
calculations for the (001) monolayer (square lattice), as- 
suming the unidirectional structure (two-fold axis of the 
molecules pointing in the direction of the external field). 
The basis for the tight-binding Hamiltonian and the hop- 
ping matrix elements were taken from the parametriza- 
tion given in Ref. 0. The splittings shown in figure 
were then used to derive an on-site coupling between the 
different orbitals. In the case of the LUMO, the on-site 
coupling is diagonal and reduces to orbital dependent 
on-site energies. The DOS for the LUMO (T lu ) and the 
HOMO (H u ) bands calculated with this model for differ- 
ent charging are shown in figure [!|] and 1101 respectively. 
We find that already for an induced charge of one carrier 
per molecule the change with respect to the unperturbed 
DOS is sizable. For q = —2 one of the Ti„-bands is al- 
ready completely separated from the other two. Also the 
HOMO density of states shows for q = 2 hardly any re- 
semblance to the original DOS, and for q = 3 also the 
-ff u -bands fall into two groups. We thus conclude that 
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FIG. 7: Splitting of the molecular levels of a Ceo molecule 
in the self-consistent multipole potential (Z < 2) for a (001) 
monolayer (square lattice) in a homogeneous external field as 
a function of the induced charge. The molecule is oriented 
such that one of its two-fold axis points in the direction of 
the external field (perpendicular to the monolayer). The top 
panel shows the positions of the split (from bottom to top) H g , 
G g , H u , Ti u , and Ti s levels. The bottom gives the splitting 
of the H u (HOMO) and T lu (LUMO) levels relative to their 
respective center of gravity. 



beyond filling |g| = 2 the electronic structure is distorted 
so much compared to the unperturbed monolayer that 
one can no longer speak of doping. 

The calculation of the density of states in a minimal 
tight-binding basis involves, of course, approximations: 
First, in the lattice, not only hopping between LUMO 
(HOMO) levels is allowed, but also hopping via energet- 
ically close-by levels. An orbital at Ae away will give a 
contribution to the hopping of about t 2 /Ae, where t is 
the hopping matrix element from the orbital, that we con- 
sider explicitly, to the orbital at Ae. The influence of this 
effect on the hopping between molecules was studied in 
Ref. and changes of the order of 5 % were found. More 
importantly, due to the deformation of the molecular or- 
bitals in the field, the hopping matrix elements between 
the HOMO or LUMO orbitals will change. For a simple 
estimate, we have performed tight-binding calculations 
of a Ceo molecule in an external homogeneous field and 
determined the average hopping matrix element between 
the tiu-orbitals following the approach of Ref. We 
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FIG. 8: Splitting of the molecular levels of a Ceo molecule 
in the self-consistent multipole potential (Z < 2) for a (111) 
monolayer (triangular lattice) in a homogeneous external field 
as a function of the induced charge. The molecule is oriented 
such that one of its three-fold axis points in the direction of 
the external field (perpendicular to the monolayer). Other- 
wise the plots are as in figure The thick lines in the lower 
panel indicate the two- fold degenerate levels. 



find that the change in the hopping matrix elements de- 
pends strongly on the orientation of the molecules. Typ- 
ical changes are of the order of 10-20%. 



VII. CONCLUSIONS 

We have analyzed the changes in the electronic struc- 
ture of a Ceo monolayer in which charge carriers are in- 
duced by the application of an external homogeneous 
electric field. We find that the effective field seen by 
each molecule in the monolayer is strongly screened, but 
that there are additional higher multipole potentials. Al- 
though these components are considerably weaker than 
the screened homogeneous field, for even I, they give a 
significant contribution to the level splitting as they are 
of first order. In addition the I = 1 (homogeneous field) 
and I = 2 potentials split the HOMO and LUMO in al- 
most the same way, so the splittings they produce add up 
or counteract, depending on the sign of the external field. 
For both, the HOMO and the LUMO, the signs turn out 
to be such, that the splitting is enhanced when the charge 
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FIG. 9: DOS (per molecule and for both spins) of the LUMO- 
band taking into account the level splitting for doping q = 
0,-1,-2. The Fermi energy is indicated. 
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FIG. 10: As in figure^lbut for the DOS of the HOMO-band 
with q = 0, 1,2,3 . 



is induced in the respective level. Thus the level that car- 
ries the field-induced charge is strongly changed by the 
effective field — a particularly unfavorable situation if 
one wants to achieve doping, i.e., filling of a level with- 
out substantially changing its electronic structure. 

There are, of course, some effects that have been ne- 
glected in our analysis: The polarizability of Cgo in- 
creases when charge is put in the T\ u orbital. Thus 
when filling the LUMO the screening of the external field 
should become somewhat more efficient. Second, due 
to the electric field the molecular orbitals are deformed 
and thus the hopping matrix elements between neighbor- 
ing molecules change. This effect depends sensitively on 



the orientation of the molecules but it typically leads to 
a slight decrease of the band width with applied field, 
making the Stark splitting even more important. Fur- 
thermore we have neglected the effect of elcctron-phonon 
coupling, which also tends to narrow the bands. This 
should be particularly important for the ff„-bandi2Si Fi- 
nally, we have not considered correlation effects. It is 
known that, at integer filling, the alkali-doped Fullerenes 
are close to a Mott transition and that they are metallic 
because of orbital degeneracy^ Due to the reduced co- 
ordination in two dimensions, a doped monolayer of Cgn 
should be even more strongly correlated, even if the or- 
bitals are still degenerate. Lifting the degeneracy, e.g., by 
the Stark effect significantly increases correlations even 
furtheriSi It therefore seems inevitable that for an inte- 
ger number of induced charges a Cgo monolayer will be 
a Mott insulator, even though it is not clear what effect 
the gate oxide might have on the Coulomb repulsion U 
between two electrons on a Cgo molecule^ 

To summarize, it seems safe to conclude that the 
physics of field-effect devices based on Ceo as active mate- 
rial should be quite different from that of the alkali-doped 
fullerenes, at least at doping levels beyond two carriers 
per molecule. 
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APPENDIX A: RESPONSE FOR GENERAL 
ORIENTATION 

To illustrate the use of the response coefficients given 
in table [HJ we show, how to calculate the response for 
a multipole potential with (7m) — (30) for different ori- 
entations of the z-axis with respect to the Ceo molecule. 
This is the first non-trivial case, as the response for mul- 
tipoles with I < 2 is isotropic. To start with, we need the 
basis functions Yi x k spanning the irreducible representa- 
tion of the Ih, that were introduced in section ILLTI For 
specific orientations, these can be found, e.g., in the Ref. 
Ill, chapter 16, or Ref.ll3l table 4.2. For arbitrary orien- 
tations, they have to be derived by explicitly finding the 
basis functions that span the irreducible representations 
of the icosahedral group. For the sake of the example, 
we consider the response for z along the 5-fold and the 
3-fold axis of the Ceo molecule. The corresponding basis 
functions are then both found in Ref. Il2l Since we are 
interested in the response to an external multipole with 
(Zm) = (30), we have to identify those basis functions 
that contain Y 30 . They are shown in table IVTl It turns 
out that for z parallel to the 5-fold axis the (Im) = (30) 
potential corresponds to a pure Ti u potential. Thus the 
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IR of I h 


5- fold axis 


3- fold axis 


Gu 




2^2 v" 1 V 
—3- 1 30 — 3 ^3c3 


T2u 


-Y 30 


— 1^30 — ^^-Yzci 



TABLE VI: Transformed I = 3 spherical harmonics, which 
are partner functions for the IR of Ih and which contain Y30. 
The second (third) column is for the case where the 3-fold 
(5-fold) axis of the molecule is parallel to the z-axis. The real 
spherical harmonic is l3 C 3 = \/2Re(V33). 



response is given by A<3 30 = a 3 3(r 2u ) V 30 . For z along 



the 3-fold axis the situation is more complicated, as the 
potential is now a mixture of partner functions of G u 
and T 2u : As can be seen from table lVll it mixes with the 
^3c3 = V / 2Re(Y33) component. By construction, the re- 
sponse matrix a in the subspace spanned by {Y30, ^303} is 
diagonal, with diagonal elements a 33 (H u ) and a 33 (T 2 u), 
when written in the basis functions listed in table IVII 
However, to obtain the multipole response we have to 
use the basis (Y 30 ,Y 3c3 ), for which there are off-diagonal 
elements: 



_ 1 f 8a 33 (H u )+a 33 (T 2u ) V&{a 33 (T 2u ) - a 33 (H u )} \ 
[033(7211) — OL3s{Hu)\ ot 33 {H u ) + 8a33(T2 ll ) / 

We thus find 

AQ30 = l[8a 33 (H u )+a 33 (T 2u )}V 30 , 
AQ 3 c3 = [a 33 (T 2u ) - a 33 (H u )] V 30 . 



APPENDIX B: COUPLING MATRICES 

In this section we discuss the calculation of the level splitting for arbitrary directions within perturbation theory 
using the coupling constants of tabic IIIII As discussed above, we restrict the analysis to I — 1 and I — 2 external 
potentials, which corresponds to x = T\ u and x = H g potentials in the icosahedral symmetry 1%. In first order 
perturbation theory, the splitting of the levels in the degenerate subspace £ nx is given by the eigenvalues of the matrix 



Hi% 2 (nx) = (nxk 2 \V eS \nxk 1 ). (Bl) 
This matrix vanishes in the case of an odd potential and the splitting is given by the second order expression 

H W {nx) = (n*k l \V cS \n'x'k')(n'x'k'\V° B \nxk 1 ) (B2) 

(n'a:') fc' 

The matrix-elements in (|B1|) and (|B2|) are given in (|llfl and involve the icosahedral Clebsch-Gordan coefficients 
Cfc 2kl (A; x 2 x\\ x). In order for the fc-indices to be defined we consider the molecule oriented with the 5-fold axis 
parallel to the z-axis (see figure This allows to label the states within a multiplet unambiguously with its C5 
index k. The ordered basis of a T lu subspace has A;-indices (0, 1,-1) whereas the ordered basis of an H u subspace is 
given by (0, 1,-1,2, —2). Note that in the case of applied I = 1 or / = 2 potential, the k index corresponds to the 
m index of the spherical harmonics. For a detailed discussion we refer to Ref. We will present the coefficients 
C% 2ki (\; X2X1; x) as matrices with respect to the indices k\ and k 2 - In order to reduce the number of matrices, we will 
give the coupling matrices for (Im) = (10), (20) potentials, which are rotated around the y-axis by an angle 9. The 
resulting matrices are then given by 

C e (X ; x 2 xi,x) = ^2 ^dfc(^) C fc (A; x 2 Xx)x), (B3) 

k 

where R l k i k (9) is the rotation matrix of the spherical harmonics in a given Z-subspace. Using the previous relations 
and equation (|11|) . the coupling matrix for an even / = 2 potential is given by 

H^(nx,9) = V 20 tx(nxnx- 1 H g )C e x {xx;Hg). (B4) 
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Note that the multiplicity label A is only relevant for the HOMO H u because H„ occurs twice in the product H U ®H U 



la 



T 2 „ © 2G g © 2H g . The coupling matrix for the odd I = 1 potential is 



H {2 \nx,9) = V 1 



10 X! 



t(nx nx' \T lu ) 2 g 



E n '. 



C (xx'\ Ti„) C (xx'\ Ti„). 



(B5) 



In the following we restrict the sum over subspaces £ n > x ' closest in energy to £ nx which is the T\ g subspace in the case 
of the LUMO and the H g and G g subspaces in the case of the the HOMO (see Fig EJ - Below, the coupling matrices 
C^(xiX2',x) which are needed to calculate the splitting of the HOMO and LUMO are given. They are traceless 
TrC e = and normalized such that Tr (C e ) T C e = 1. The coupling matrices which describe the splitting of the LUMO 
are 



C (Ti u Tig;Ti z 



C e (Ti u Ti u ; H g ) — — 




— sin 9 — sin 9 
-V2cos6> 
\/2cos( 



V 



4 , - - 2 V6 sin 2 
\/3sin26> 
-\/3sin26l 



V3sin26» 

\/6 sin 2 
a/6 sin 2 



-\/3sin20 \ 
V6sin 2 6> 



(B6) 



(B7) 



The eigenvalues of these two matrices are independent of 9 and given by (0, — , ) and (— - 2 =, -i=, -i=) which 
implies that the splitting is independent of the orientation of the molecule with respect to the direction of the applied 
/ = 1 and / = 2 potentials. The coupling of the HOMO (H u ) to the lower lying H g and G g levels is given by 



C e {H u H g ; T\ u ) 



C e (H u G g ;Ti u ) 



2\/5 



2x75 



/ 





a/3 sin 


— a/3 sin 







-a/2 cos 





a/2 sin^ 




— a/3 sin 





a/2 cos 










a/2 sin 





-2a/2cos 


V 








-\/2sm0 







\/3sin0 


2a/2cos0 





-71 sin 




—a/3 sin 





-2a/2cos0 










-\/2sin0 





—a/2 cos 


V 








a/2 sin 9 


-2= sing 

\/2 






-A/2sim 


2a/2~COS( 







v/2 
3 



smt 



a/2cos0 / 



(B8) 



(B9) 



The coupling of the HOMO among themselves is given by 



C (1; H U H U ; H g ) 



C (2;H u H u ;H g ) 



4a/5 



12 sin^ 
f sin 29 
| sin 29 



| sin 2 9 



sin 



2 



sin 29 
| sin 26» 



| sin 2 f 
| sin 2 f 



n 29 

2-3sin 2 -3 sin 2 9 
-3 sin 2 (9 2-3sin 2 
-3 sin 20 
3 sin 2 9 



sin 29 \/^sm~ 
-3 sin 20 
3 sin 2 9 
3 sin 2 9 2 - 3 sin 2 9 
-3 sin 20 3 sin 20 



f sin 29 



| sin 2 f 



■ sin 29 

2-3 sin 2 9 sin 2 9 sin 29 

sin 2 9 2 - 3 sin 2 9 sin 2 9 

sin 20 sin 2 -2 + 3 sin 2 

sin 2 9 sin 29 sin 29 



sin o 
3 sin 2 9 
-3 sin 29 
3 sin 20 
2-3 sin 2 6 



| sin 2 i 



sin 2 (9 
sin 29 
sin 29 
-2 + 3 sin 2 9 



(BIO) 



(Bll) 
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In equation (|B5H the product C 9 (xx'; Ti u ) T C 9 (xx'; T\ u ) enters. For the matrices given in (|B6(1 . (|B8|I and i|B9|) these 
products can be expressed in terms of the I = 2 coupling matrices l|B7|l . (|B10|I and (|B11() : 



C e (Tx u Ti g \ Ti u ) T C e (Ti u Ti g ; T\ u ) — — + — j= C 9 (Ti u Ti u : H g ), 



C e (H u H g ;Ti u ) T C (H u H g ]Ti u ) 
C e (H u G g ; Ti u ) T C 9 (H u G g ; T\ u ) 



- + C e (l; H U H U \ H g ) — — C 9 {2; H u H u ;H g ), 

- + C e (l; H U H U ; H g ) + — C 9 {2; H u H u ;H g ). 



(B12) 

(B13) 
(B14) 



Using these relations, the total coupling matrix (neglecting the constant terms in the previous relations) due to the 
applied I = 1 and I — 2 potential, is given by 



H'(T lu ,6) = [V^c 1 + V w c 2 ) C e (T lu T lu ;H g ), 



H'(H U ,6) = V? di cos Si C 9 (1; H u H u ;H g )+ sin S^ 9 (2; H u H u ;H g ) 



V 20 d 2 



cos 5 2 C 9 {1; H U H U - H g ) + sin 5 2 C 9 (2; H U H U ; H g ) 



(B15) 
(B16) 
(B17) 



where ci = -4= 



_ 1 t(T lu T lg ;T^ 



^ — _ Et — and c 2 — t(Ti u Ti u ; H g ). Similarly we have d\ cos5\ = 



t(H u H g \T lu 



t(H^G g ;T lu ) 2 
E Hu -E Gg 



d% sin S% = 



t(H u H g ;T lu y t(H u G g ;T lu y 
Eh u — Eh q Eh u — Eq 



and d 2 cos<5 2 = ti(H u H u ; H g ), d 2 sin<5 2 = t 2 (H u H u ; H g ). Equation (|B15() 
implies that the contributions of the / = 1 and I = 2 to the splitting of the LUMO add up trivially. Using the values 
in table ITTT1 and the energies of figure EH yields the values 8\ = 0.064 and S 2 = 0.037, which are almost equal when 
compared to 7r. This can be understood by the following remarks: <5i — in the case of t(H u H g ; T\ u ) — t(H u G g ; T\ u ) 
and En g = Ec g - Furthermore, it can be shown that 5 2 — arctan(l/Vl25) « 0.090 assuming that the angular 
dependence of the HOMO is given by I = 5 spherical harmonics. Taking an average value of 5 — 0.050 yields the 
approximate relation 



H'(H U ,6) a (V* di + V 20 d 2 ) [cos 6 C 9 (I; H u H u ;H g ) + smSC 9 (2; H U H U ; H g ) 



(B18) 



which shows that, to a good approximation, the contributions of the I = 1 and I = 2 potential to the splitting of the 
HOMO add up trivially. 



APPENDIX C: CALCULATION OF THE MATRIX f3 

In this section it is shown how to calculate the matrix (3 appearing in equation l|l(jf> . The second term on the 
right side of this equation describes the coefficients of the term X^r^o ^'( r — R) w hich enters the screened potential 
V~ scr and which describes the potential induced by all neighboring sites. Using the definition Q), we can rewrite this 
expression as 

£ V\V ~Ri) = J2 Qlrn £ IL(T ~ Ri) = Yl ® lm E {-VTIl-n (' - Hi)- (ci) 
R;#0 Im Ri^O Im R;#0 

The function I; m (ri — r 2 ) can be decomposed using the translation formula (for r± < r 2 ) 

°° / (21 4- 1)1 

Wri-r 2 )= (~ 1 )' 2 V r 2 LiW2*V S ^m^^^m^ri)^^^), (C2) 

h,'2 = « V ^ ' )-\ U' mijm2 

! 2 -( 1 =i 

where C, L ^f , ^ denote the Clebsch-Gordan coefficients. This formula can be found in different forms in the literature, 
see for example Ref. Il4ll5l Substituting l)C2jl in the sum (|C1|) with ri = r and r 2 = Ri yields 

U'(r-Ri) = 2J PhmthT^Qhn^Rhnntr), (C3) 
Ri^O ! i"«2 

!2 m 2 
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where the matrix (3 is given by 



Phm 1 l 2 m2 ^ I X J y (2^i)!(2Z2)! ™ 2 / ~ ^ J /i+i2 "11^2 l-^iJ 



f _,™ !+ i a / {l 1 + l 2 +m 1 -m 2 )\(l 1 +l 2 -m 1 + m 2 )\ ^ 

In the last equality in (|C4ll the explicit form of the Clebsch-Gordan coefficients was usedii One verifies, that the 
matrix (3 is complex conjugate under the exchange of all indices. The remaining sums over the lattice sites in (|C4|) 
are easily performed by computer. 
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